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ABSTRACT 

We introduce a simple model to explore the star formation histories of disk galaxies. We 
assume that the disk origins and grows by continuous gas infall. The gas infall rate is parametrized 
by the Gaussian formula with one free parameter: infall-peak time t p . The Kcnnicutt star 
formation law is adopted to describe how much cold gas turns into stars. The gas outflow process 
is also considered in our model. We find that, at given galactic stellar mass M„, model adopting 
late infall-peak time t p results in blue colors, low metallicity, high specific star formation rate 
and high gas fraction, while gas outflow rate mainly influences the gas-phase metallicity and 
star formation efficiency mainly influences the gas fraction. Motivated by the local observed 
scaling relations, we construct a mass-dependent model by assuming low mass galaxy has later 
infall-peak time t p and larger gas outflow rate than massive systems. It is shown that this model 
can be in agreement with not only the local observations, but also the observed correlations 
between specific star formation rate and galactic stellar mass SFR/M* ~ M* at intermediate 
rcdshift z < 1. Comparison between the Gaussian- infall model and exponential-infall model is 
also presented. It shows that the exponential-infall model predicts higher star formation rate at 
early stage and lower star formation rate later than that of Gaussian-infall. Our results suggest 
that the Gaussian infall rate may be more reasonable to describe the gas cooling process than 
the exponential infall rate, especially for low-mass systems. 

Subject headings: Galaxies:evolution - Galaxies:photometry - Galaxies: stellar content 



1. Introduction 

More than thirty years ago, Visvanathan & Gri- 
ersmith (1977) found that early-type spiral galax- 
ies (SO/a to Sab) exhibit optical color-magnitude 
relation (CMR), that is, luminous spirals show 
redder colors than less-luminous systems. 

Later, many studies confirmed this phenom- 
ena and extended the CMRs of spirals from op- 
tical to infrared bands (Tully ct al. 1982; dc 
Jong 1996; Mobasher et al. 1986; Tully et al. 
1998; Peletier & de Grijs 1998; Bell & de Jong 
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2000; Blanton et al. 2003; Baldry et al. 2004; 
Macarthur et al. 2004). Chang et al. (2006a) col- 
lected a sample of late-type galaxies from Sloan 
Digital Sky Survey (SDSS) and Two-Micron All- 
Sky Survey (2MASS) and found that, after cor- 
recting the effect of dust attenuation, the optical- 
infrared (optical-IR) colors present tighter corre- 
lations with the absolute magnitude than optical 
colors. Based on a matched sample of Galaxy 
Evolution Explorer (GALEX) and SDSS, Wyder 
et al. (2007) and Schiminovich et al (2007) ex- 
plored the ultraviolet-optical (UV-optical) CMRs 
for both disk and bulge-dominated galaxies. They 
also found that disk-dominated galaxies show a 
tight correlation between specific star formation 
rate (SFR/M*) and galactic stellar mass (M»), 
which is denoted as the mass-SFR relation for con- 



venience. 
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On the other hand, interpretation of these ob- 
served CMRs of disk galaxies is still an open ques- 
tion, van den Bosch (2002) and Bell ct al. (2003) 
found that, contrary to the observed trends, their 
semi-analytic models (SAMs) predicted that faint 
spirals should be slightly redder than bright spi- 
rals. By introducing additional mechanics like ac- 
tive galactic nucleus (AGN) feedback, the more 
recent SAMs can predict positive slopes of the 
CMRs of spirals, i.e., lower mass systems being 
bluer than massive systems (e.g. Kang et al. 2005, 
Croton et al. 2006). However, the slopes of the 
predicted CMRs are still shallower than that of the 
observations (Kang et al. 2005). This is mainly 
due to the fact that it is difficult for the SAMs 
to build a scenario that permits low mass spiral 
to be averagely younger than high mass system in 
the hierarchical cosmogonies. 

Parametrized modelling, which adopts phe- 
nomenological descriptions for some complicated 
processes, is another fruitful tool to explore the 
galactic formation and evolution (Tinslcy 1980). 
In fact, parametric models have been successfully 
applied to investigate the formation and evolu- 
tion of the Milky Way disk (Chang ct al. 1999; 
Chiappini ct al. 2001; Rcnda ct al. 2005; Naab 
& Ostrikcr 2006; Fu ct al. 2009; Matteucci et 
al., 2009) and other nearby disk galaxies (Prant- 
zos& Boissier 2000; Molla & Diaz 2005; Dalcan- 
ton 2007; Yin et al. 2009). Besides these spe- 
cific galaxies, Boissier & Prantzos (2000, hereafter 
BP00) extended their calibrated model on the 
Milky Way disk to study the chemical and spectro- 
photometric evolution of disk galaxies and com- 
pared their model predictions with the observed 
CMRs of de Jong (1996). They found satisfactory 
agreement with local observations provided that 
the star formation time-scales of low-mass galax- 
ies are longer than their high-mass counterparts. 
However, their model predictions have not been 
compared with observations at higher rcdshift due 
to the limited observations at that time. 

Recently, several new scaling relations of disk 
galaxies are statistically found and huge amount 
of observed data are available, not only in the local 
universe but also for intermediate and high red- 
shift galaxies. Motivated by the work of BP00, we 
try to construct a new parametric model to ex- 
plain the observed correlations in both local and 
intermediate redshift universe. Our results show 



that the mass-SFR relations at intermediate red- 
shifts provide tight constraints on the modelling of 
gas cooling history. Our results suggest that the 
Gaussian infall rate may be more reasonable to de- 
scribe the gas cooling process than the exponential 
infall rate, especially for low-mass systems. 

The outline of this paper is as follows. In Sec- 
tion 2, we describe main assumptions and ingredi- 
ents of our model in details. Comparisons between 
model predictions and observations and our main 
results are given in Section 3. Discussions are pre- 
sented in Section 4 and the last section summarizes 
our main conclusions. 

2. The model 

We assume that the disk galaxy has been em- 
bedded in a dark matter halo and primordial gas 
within the dark halo cools down to form a rota- 
tionally supported disk. The disk is basically as- 
sumed to be sheet-like and composed by a set of 
independent rings, each 500pc wide, that is, there 
is no gas radial flow between different rings. The 
details and essentials of our model arc outlined in 
below sections. 

2.1. Present-day disk galaxy 

In the local Universe, galactic stellar mass M* 
is the primary parameter that determines main 
properties of the galaxy (Kauffmann et al. 2003a, 
2003b), including the colors of disk galaxy (Chang 
et al. 2006a). In this paper, we assume M* is the 
primary parameter and all free parameters in our 
model are dependent on , where M* is referred 
to the stellar mass of a galaxy at present time (i.e., 
z=0) if there is no specific notations. 

For a present-day galaxy with galactic stellar 
mass M* , we assume the stellar mass surface den- 
sity in the current time £*(r, t g ) follows an expo- 
nential profile: 

£*M 9 ) = £.(0,< g )exp(-r/r d ), (1) 

where t g is the cosmic age at the present time, 
which is 13.5 Gyr for the cosmogony adopted here 
(see following subsection). E^Ojig) is the central 
stellar mass surface density in the present-day and 
rd is the disc radial scale-length. Since a pure 
exponential disk has: 

M* = £* {0,t g )2mi, (2) 
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£*(r, t g ) can be obtained if M„ and the M* ~ 
relation are given. 

Shen et al. (2003) collected a large sample of 
SDSS galaxies and obtained the mean M* ~ R$o 
relations for both early-type and late-type galax- 
ies in the local universe, where -R50 is the aperture 
that contains 50% of Petrosian flux (Stoughton et 
al. 2002). We adopt their observed M* ~ i?so 
relation of late-type galaxies in z-band (equation 
18 and table 1 of Shen et al. 2003) to determine 
Td{M^) by adopting = i?so/1.68 for a pure ex- 
ponential disk, since z-band is a good indicator of 
the stellar mass. 

2.2. Gas infall rate 

Gas cooling is one of the key processes in mod- 
eling galactic formation and evolution. Semi- 
analytical models, which have made large suc- 
cesses in explaining a lot of observations (Kauff- 
mann et al. 1993; Somerville & Primack 1999; 
Cole et al. 2000; Kang et al. 2005; Croton et 
al. 2006; Fontanot et al. 2009), have presented a 
method to describe gas cooling process. However, 
there are still some open questions, for example, 
how to describe the mass and angular momentum 
distributions of baryons in dark matter halo and 
their variations with time, and whether the spe- 
cific angular momentum conserves or not during 
gas cooling processes etc. (Efstathiou 2000; Bul- 
lock et al. 2002; Chen et al. 2003). None of these 
questions has an answer clear enough to construct 
a simple model, at least at the present stage of our 
knowledge. 

In this paper, we turn to another approach, i.e., 
the phcnomcnological approach, which adopts a 
parametric description for some complicated pro- 
cesses. Comparisons between model predictions 
and observations may suggest possible ranges for 
model free parameters in the model and then help 
to know more information about main properties 
of the processes described by these parameters. 
Wc assume that there is 1 Gyr time delay for the 
disk galaxy formation, which corresponds to z ~ 6 
under the adopted cosmology in this paper. After 
that, the disk origins and grows by continuous gas 
infall from the dark halo and the metallicity of the 
infalling gas is primordial. Differing from BP00, 
who adopted an exponential infall rate, we follow 
Prantzos & Silk (1998) and Chang et al. (1999) 
and adopt a Gaussian formula of gas infall rate. 



Comparison between exponcntial-infall model and 
Gaussian-infall model will be presented in Section 
4. 

The Gaussian infall rate /m(r, t) (in units of 
M Q pc~ 2 Gyr~ 1 ) is assumed to be: 

f in{r ,t) = ^fle-^ 2 /^ (3) 

where t p is the characteristic time corresponding 
to the maximum gas infall rate (hereafter, wc call 
t p as the infall-peak time) and a is the full width 
at half maximum of the peak. We adopt a = SGyr 
throughout this paper to promise a relatively wide 
distribution of the infall rate. In fact, our calcula- 
tions show that the varying of a in the range 2^4 
Gyr do not have large influence on the results. The 
A(r) are a set of separate quantities determined by 
the boundary condition £*(r, t g ). In practice, we 
itcratively obtain A(r) by normalizing the resulted 
stellar mass surface density distribution at t = t g 
to be very close to its observed value £*(r, t g ). 

We should emphasize that the infall-peak time 
t p is the most important free parameter in our 
model that determines main properties of gas in- 
fall history. Indeed, t p regulates the shape of gas 
accretion history. For example, t p — > corre- 
sponds to a time-declining infall rate, while t p — > 
00 corresponds to a time-increasing gas infall rate. 
Since we only focus on the global quantities of 
galaxy in this paper, wc assume t p does not vary 
with radius r but may be dependent on the galac- 
tic stellar mass M*. 

In this paper, we assume a cosmology with 
n M = 0.3, 0\ = 0.7 and Hq = TOkms^Mpc" 1 . 
Correspondingly, the cosmic time at z =0 is t g = 
13.5 Gyr and the cosmic time at z =3 (0.4) is t = 
2Gyr (9Gyr). 

2.3. Star formation law and stellar evolu- 
tion 

Star formation (SF) process is so complicated 
that it is very difficult to fully understand the 
physical nature of star formation regulation (Silk 
1997, Kennicutt 1998a,b, Elmegreen 2002). For- 
tunately, observations of galaxies on global scales 
reveal a tight correlation between the average star 
formation rate per unit area and the mean surface 
density of cold gas (e.g. Kennicutt 1998a, b). This 
observed relation is later called as Kennicutt SF 
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law. Recently, Kennicutt et al. (2007) found that, 
on spatial scales extending down to at least 500 
pc, the star formation rate (SFR) surface density 
also correlates with the local gas surface density, 
following a similar power law as that of the global 
quantities. 

We adopt the Kennicutt SF law in this paper: 

f(r,t) = 0.25S s (r,t) 1 - 4 , (4) 

where the SFR surface density ^(r, t) is measured 
in units of M Q pc _2 Gyr _1 and the cold gas surface 
density T, g (r,t) is in units of M Q pc~ 2 . 

Regarding the chemical evolution of the galactic 
disk, both instantaneous-recycling approximation 
(IRA) and instantaneous mixing of the gas with 
stellar ejecta is assumed, i.e., the gas in a fixed ring 
is characterized by a unique composition at each 
epoch of time. We take the classical set of equa- 
tions of galactic chemical evolution from Tinsley 
(1980) and adopt the return fraction R = 0.3 and 
stellar yield y = IZq, 

The updated stellar population synthesis (SPS) 
model of Bruzual & Chariot (2003), i.e., CB07, is 
adopted in the present paper, with the stellar IMF 
being taken from Chabrier (2003). The lower and 
upper mass limits are adopted to be O.IMq and 
IOOMq, respectively 

2.4. Gas outflow rate 

Gas outflow process may also influence the evo- 
lution of galactic disk. During the continuous 
star formation and evolution process, supernova 
(SN) explosions eject part of energy into interstel- 
lar medium and some amount of gas will be heated 
and expelled from the galactic disk if its kinematic 
energy exceeds the bonding energy. Therefore, gas 
outflow rate may be proportional to the SN explo- 
sion rate and then proportional to SFR surface 
density $(r, t) if the initial mass function (IMF) 
does not vary with time. Furthermore, gas out- 
flow rate may also correlate with galactic mass 
since low mass galaxies have shallower potential 
wells and then it is much easier for them to expel 
gas out than high mass systems (Kauffmann et al. 
1993). 

In this paper, we assume the outflowing gas 
has the same metallicity as that of the interstellar 
medium at that time and will not fall again to the 
disk. The gas outflow rate f ou %{r,t) (in units of 



M Q pc 2 Gyr ) is assumed to be proportional to 
the SFR surface density *(r, t): 

/out(M) = b out $>(r,t), (5) 

where 6 out is another free parameter in our model 
and may also be a function of galactic stellar 
mass M*. Integrating the left side of the above 
equation, we derive J Q 9 / ou t(r, t) drdt = M out , 
where M out is the total mass expelled from the 
galaxy. Integrating the right side of Eq. 5 gets 
So' So°° & out*(r, t) drdt = 6 out (l - R)M„ where R 
is the return fraction. Then, we derive M out = 
&out(l — R)M*. Therefore, 6 ou t is proportional to 
the fraction of total mass expelled relative to the 
final mass locked in stars. 

In summary, our primary input parameter is 
the galactic stellar mass M*. Free parameters are 
infall-peak time t p and outflow parameter b ou t- 
The combination of gas infall rate and outflow 
rate determines the variation rate of the total 
(gas+star) mass surface density. Assuming the 
initial total surface density to be zero, we can 
numerically obtain the total mass surface den- 
sity at any time after free parameters being given. 
Adding the SF law, we can describe how many 
cold gas turns into stellar mass and then calcu- 
late the chemical and color evolution of the disk. 
For simplicity, we assume that there is no scat- 
ter of the model parameters t p and 6 ou t at any 
given M*. In other words, our simple model only 
investigates the average trend of the correlations 
between different quantities and does not include 
their scatters. 

3. Model results versus observations 

Using the above model, for one galaxy at given 
M*, we independently calculate the evolution of 
cold gas in each ring along the galactic disk. At 
each time step, we derive a set of individual quan- 
tities, such as the SFR surface density, gas mass 
surface density, multi-band luminosity and metal- 
licity, and then integrate them to obtain the global 
quantities for the whole galaxy. After that, our 
model predictions can be directly compared with 
observations. 

In this section, we present our main results step 
by step. Firstly, we briefly describe the adopted 
observations to be compared with model predic- 
tions. Then, we investigate the influence of free 
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parameters on model predictions, including not 
only t p and 6 ou t, but also the Constance in the 
SF law. Finally, we construct a reasonable mass- 
dependent model and explore if its predictions can 
be consistent with the observed data. 

3.1. Observations 

Here, we show and summarize the observed 
mass-color relations, the mass-mctallicity relation, 
mass-SFR relation and the mass-gas fraction re- 
lation for local galaxies, which are mainly taken 
from papers based on SDSS database. We also in- 
clude the observed mass-SFR relations at redshift 
z = 0.3,0.5,0.7,0.9. 

Mass-color relations: Based on the matched 
sample of SDSS and GALEX, Salim ct al. (2007) 
and Schiminovich et al. (2007) obtained the UV- 
optical color magnitude diagrams of star forming 
galaxies in the local universe and investigated their 
physical properties as a function of specific star 
formation rate {SFR/M*) and stellar mass (M*). 
We take the observed (NUV — r) ~ M* relations 
of disk-dominated galaxies from Table 2 of Schimi- 
novich et al. (2007) and plot them in the left- 
upper panel of Fig. 1, where the open circles are 
the data without dust-correction and the squares 
denote those after dust corrections. The error-bars 
denote the dispersions of the data distributions. 

Chang et al. (2006a) collected a large sample of 
late-type galaxies from SDSS (DR4) and 2MASS 
to explore their CMRs in both optical and optical- 
IR bands. The photometric aperture mismatch 
between SDSS and 2MASS is taken into account 
by correcting the SDSS magnitudes to the isopho- 
tal circular radius, where the 2MASS magnitudes 
are measured, based on the radial profile of surface 
brightness of SDSS galaxies. We take the observed 
multi-band (optical and optical-IR) color- M* rela- 
tions from Chang et al. (2006a) and plot them in 
Fig. 1, where the square denotes the mean value 
of color after dust-corrections and the error-bar 
represents the distribution dispersion in each bin. 
Since the SPS models seem to be not accurately 
reproducing the shape of the spectrum between 
4000 and 5000 A, which causes an offset in g-r and 
r-i as the r filter passes through this (rest-frame) 
wavelength region (Wake et al. 2006), we ignore 
the r filter data and replace them with g filter data 
in this paper. 



Mass-metallicity relation: The mass-metallicity 
relation is another tight constraint on models of 
galactic formation and evolution. Based on the 
SDSS imaging and spectroscopy data, Tremonti 
et al. (2004) derived gas-phase oxygen abun- 
dances and stellar masses of a large sample of 
star-forming galaxies and found a tight correla- 
tion between stellar mass and metallicty, span- 
ning over 3 orders in stellar mass and a factor of 
10 in metallicity. Liang et al. (2007) adopted an- 
other method to measure the gas-phase metallicity 
of a similar sample of star-forming galaxies from 
SDSS (DR4). Within 27 stellar mass bins from 
log^M^/Mo) = 8.0 to 10.6, they stacked the spec- 
tra of several hundreds (even thousands) galaxies 
in each bin and derived their electron temperature 
T e and hence obtained the T e -based oxygen abun- 
dance. In this paper, we adopt these two sets of 
observed relations between gas phase metallicity 
log(0/H) + 12 and stellar mass Af». The results 
are plotted in the right-second panel of Fig. 1, 
where the data taken from Tremonti et al. (2004) 
are shown by yellow shaded area and that taken 
from Liang et al. (2007) are denoted by green 
filled circles with error bars. 

Mass-SFR relation: The specific star for- 
mation rate SFR/M* is an important quantity 
describing the star formation history. If the ini- 
tial mass function does not vary with time, stel- 
lar mass should be roughly proportional to 
mean star formation rate during the whole history 
< SFR >, thus, SFR/M* cx SFR/ < SFR >. 
Therefore, a galaxy having higher SFR/M* means 
larger fraction of its stars formed in later time. 
The observed SFR/M* ~ Af* relations of star- 
forming galaxies in the local universe are taken 
from Schiminovich ct al. (2007) and plotted in 
right-third panel of Fig. 1, where the SFRs taken 
from Brinchmann et al. (2004) and from Salim 
et al. (2007) are denoted as the open circles and 
squares, respectively. 

Besides the local data, the observed mass-SFR 
relations at different redshift are also available. 
Zheng ct al. (2007) measured the SFR from UV 
and IR luminosities for a sample of ^15,000 galax- 
ies selected from two COMBO-17 fields. Combin- 
ing with the stellar mass of these galaxies derived 
by Borch et al. (2006), Zheng et al. (2007) derived 
the observed mass-SFR relations at four redshift 
bins z — 0.3, 0.5, 0.7, 0.9. We adopt their data and 
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Fig. 1. — Comparisons between model predictions 
and local observations. Different lines correspond 
to various parameter groups: dot lines {t pi b out ) 
= (5Gyr,0), solid lines (t p ,b out ) = (lOGyr.O), 
dot-dash lines (t p ,b out ) = (10Gyr,l). Dash lines 
adopts (tpjbout) = (10Gyr,0), but a modified star 
formation law \I/(r, t) = 0.05£g' 4 (r, t) is assumed 
to reduce the star formation efficiency. The ob- 
served data are described in Section 3.1 in details. 

plot them as filled squares in Fig. 2. 

Mass-gas fraction relation: The observed 
gas fraction is also important constraint on the 
model. Zhang ct al. (2009) collected a sample of 
721 galaxies from cross-matching the SDSS (DR4) 
galaxies, the 2MASS XSC and the HyperLeda HI 
catalogue (Paturcl ct al. 2003) and obtained their 
observed gas fraction. The observed data of gas 
fraction are taken from sample 1 of Zhang et al. 
(2009) and plotted as small dots in the right-fourth 
panel of Fig. 1. In the same panel, the filled tri- 
angles represent the observed data of gas fraction 
and stellar mass for a sample of 60 local galaxies, 
which are taken from Table 1 of McGaugh (2005). 

3.2. Influence of free parameters 

In this section, we show model predicted scal- 
ing relations with different settings of free parame- 
ters and compare them with the observations. Fig. 



Fig. 2. — Model predictions compare to the ob- 
served data at redshift z ~ 0.3, 0.5, 0.7, 0.9, re- 
spectively. The notations of different lines are the 
same as Fig.l. The observed data are taken from 
Zheng et al.(2007). 

1 compares model predictions with local observa- 
tions. Fig. 2 plots the observed and model pre- 
dicted mass-SFR relation at z ~ 0.3,0.5,0.7,0.9, 
respectively. The observed data in Fig. 1 and 

2 are described in details in subsection 3.1. We 
emphasize that, for both observations and model 
predictions, M*(i) in Fig. 2 refers to the galactic 
stellar mass at given redshift (or evolutionary time 
step). 

We start from a prior model adopting t p = 
lOGyr and b ou t = 0, which corresponds to a late 
infall-peak time and no gas outflow is considered. 
The results of our prior model arc plotted as solid 
lines in Fig.l and 2. It should be pointed out that 
the construction of this prior model is arbitrary 
and far from a reasonable description of the for- 
mation of any galactic disk. This setting only aims 
to provide a benchmark to be compared with other 
model predictions and then investigate the influ- 
ences of different settings of free parameters on 
model results. 

Firstly, we fix b out = and change t p from 
lOGyr to 5Gyr, i.e. replacing the prior model with 



a shorter infall-peak time. The model predictions 
are shown as the dot lines in Fig. 1 and 2. It can 
be seen that the model adopting later t p predicts 
bluer color, lower metallicity, higher gas fraction 
and higher specific star formation rate SFR/AI* 
both in the local Universe and at intermediate red- 
shift, at least up to z ~ 1. This is mainly due to 
the fact that, in our model, the setting of later 
infall-peak time corresponds to slower gas infall 
process and thus slower star formation processes. 
Furthermore, Fig. 1 and 2 show that the separa- 
tions between the solid and dash lines almost cover 
the observed range of different data sets, especially 
the SFR/M* at intermediate rcdshifts. In other 
words, the infall-peak time t p significantly influ- 
ence the galactic evolution. 

The impact of outflow (described by b ou t) on 
the model predictions is presented by the dot-dash 
lines in Fig. 1 and 2, where we adopt 6 ou t = 1.0 
and t p = lOGyr. Comparison between solid lines 
and dot-dash lines shows that the gas outflow pro- 
cess mainly influences the final gas-phase metal- 
licity, since it takes a fraction of metal away from 
the disk. Moreover, after including gas outflow 
process, the model predicts bluer optical-IR (both 
g — J and g — K) color, while the impact on the 
optical color is very small, if there is any. This is 
mainly due to the fact that the stellar mean age 
influences the optical color more strongly than the 
optical-IR color, while the mean stellar metallicity 
influences the optical-IR color more strongly than 
the optical color (Chang et al. 2006b). 

The star formation law is another important in- 
gredient in our model. In this paper, we adopt a 
Kcnnicutt star formation law. Although it works 
well in the region of high surface density, we should 
be caution to apply it to low density region, such 
as the outskirts of disk galaxies. To explore the 
influence of the adopted star formation law on our 
model predictions, we adopt a modified star for- 
mation law *ff(r,t) — 0.05£g 4 (r, t) to reduce the 
star formation efficiency (SFE) and plot the model 
predictions as dash lines in Fig. 1 and 2. The free 
parameter combination is adopted to be (t p ,b ou t) 
— (10Gyr,0). It shows that the low SFE results in 
bluer color, lower metallicity and higher SFR/M*. 
The most outstanding point appears in the mass- 
gas fraction relation. The model adotping the low 
SFE predicts much higher gas fraction than that 
of standard SFE. In other words, the predicted 




log(M./M ) 

Fig. 3. — The adopted infall-peak time t p (M*) 
and outflow parameter b out (M*) in our mass- 
dependent model. 

gas fraction is very sensitive to the adopted SFE. 
However, except the panel of gas fraction, the sep- 
arations between the solid and dash lines are very 
small compared to the observed range of different 
data sets, i.e., the variation of SFE may not dra- 
matically change our model predictions. Consid- 
ering the fact that there still lacks an alternative 
choice more reasonable than the Kennicutt star 
formation law, we choose to adopt the Kennicutt 
star formation law throughout this paper if there 
is no specific statement. 

3.3. A mass-dependent model 

Model calculations in above section show that 
galaxies form by later infall-peak time will result 
in bluer color and higher specific star formation 
rate, while gas outflow process will decrease the 
gas metallicity in the present day. The local obser- 
vations show that less massive disk galaxies have 
bluer color, lower metallicity when compared with 
their high mass counterparts (see Fig. 1). There- 
fore, it is naturally to assume that low mass disk 
galaxies may have later infall-peak time and larger 
outflow parameter than massive systems. 

Motivated by the above discussions, we try to 
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Fig. 4. — Our mass-dcpcndcnt model predic- 
tions (solid lines) compare to the local observa- 
tions. The dash lines show the predictions of 
the exponcntial-infall model, which is described in 
Section 4 in details. 

construct a model in which the infall-peak time 
t p and outflow parameter bout are only depen- 
dent on galactic stellar mass M*. The adopted 
t p (M*) and b out (M*) are plotted in Fig. 3. The 
adoption considers the facts below. Since mas- 
sive spirals have very old stellar population and 
their potential wells are so deep that there is al- 
most no outflow occurs, we set one boundary con- 
dition at high-mass end logM^/MQ = 12 to be 
t p = 2Gyr and 6 ou t = 0. Regarding the facts that 
low mass spirals are observed to have blue colors 
and their potential well is shallow and then should 
have large outflow rate, another boundary condi- 
tion of t p = ll.SGyr and 6 out = 0.7 is chosen at 
the low-mass end logM*/M Q = 8. 

Moreover, the general trends of t p (M*) and 
bout(M^) are consistent with previous studies 
of disk galaxies in the local group. For exam- 
ple, the Milky Way disk have stellar mass about 
logM*/M Q = 10.6 (Yin ct al. 2009). Our mass- 
dependent model adopts the corresponding infall- 
peak time about 7.8 Gyr, which is roughly in 
agreement with our previous chemical evolution 
models (see Chang et al. 1999 and Yin et al. 2009 



Fig. 5. — Comparison between mass-dcpcndcnt 
model predictions and observations at intermedi- 
ate redshifts. The notation of solid and dash lines 
are the same as Fig. 4. 

for details). For M33, chemical evolution models 
show that this low mass disk galaxy should form 
by slow accretion process (Wang et al. 2009; Ma- 
grini ct al. 2007), which equals to a late infall-peak 
time in our model . 

The mass-dependent model predictions are 
plotted as solid lines in Fig. 4 and 5. As ex- 
pected, under the assumption that t p and b ou t 
all decrease with the increase of M*, our model 
do predict mean trends of most of observations 
both in local Universe and at different redshift, 
at least to z ~ 1. It should be emphasized that, 
although the accurate values of free parameters 
in our model arc not unique for any individual 
galaxy, the main trend that both t p and b ou t de- 
crease with the increase of M* are robust. 

4. Discussions: exponential-infall model 

In above section, we adopt the Gaussian in- 
fall rate. However, the exponential infall rate 
is much more widely assumed in previous mod- 
els. In this section, firstly, we will compare our 
mass-dependent model assuming the Gaussian in- 
fall rate with that adopting the exponential infall 
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Fig. 6. — Gas infall history and total star forma- 
tion history of galaxies of logM„/MQ = 9.0 and 
logM*/M Q = 12.0. The y-axis of two left pan- 
els are normalized by 10 9 M Q Gyr~ 1 , while that of 
two right panles are normalized by 10 12 M Q Gyr _1 . 
The solid and dash lines represent the model re- 
sults adopting Gaussian and exponential infall 
rate, respectively. 

rate. We hope the observed data may provide con- 
strains on the form of gas infall rate and then help 
us to know more about the gas cooling process in 
disk galaxies. Secondly, we will discuss the differ- 
ence of our model with that of BPOO. 

In order to investigate the possible influence of 
the infall form on the galactic evolution, we fix 
other ingredients of the mass-dependent model but 
replace the Gaussian infall rate by an exponential 
one: 

g m (r,t) = B{r)e- l l\ (6) 

where gi n (r, t) is the gas infall rate in units of 
M0pc~ 2 Gyr -1 , r is the infall time-scale. The 
B{r), which is similar as A(r) in Eq. 3, is a 
set of normalization quantities constrained by the 
boundary condition £*(r, t g ). 

Accordingly, we also assume that infall time- 
scale r is only dependent on the stellar mass. To 
be consistent with the observed data, we adopt the 



form of r(M„) as: 



G yr e (logM»/M o -10) _ e -(logM»/M -lO) ' V> 

which gives that the infall time-scale of massive 
galaxy (Af» > 10 10 M Q ) is positive and decreases 
with the increase of M* , while the infall time-scale 
of low-mass galaxy (M» < 10 10 Mq) is negative 
and increases with the decrease of M*. It should 
be pointed out that, although r(M t ) is discontinue 
at M* = 1O 1O M0, the gas infall history described 
by the adopted r(M*) varies smoothly from in- 
creasing with time to being constant and then to 
decreasing with time as galaxy changes from low- 
mass to intermediate and then to high-mass sys- 
tems. 

To visually compare the Gaussian and expo- 
nential infall rate, Fig. 6 plots the gas infall 
history and total star formation history of two 
galaxies A/* = 10 9 M Q and M* = 10 12 M o , which 
are normalized by lO 9 M0Gyr _1 (left two pan- 
els) and lO 12 M0Gyr _1 (right two panels), respec- 
tively. The solid and dash lines represent the 
model results adopting Gaussian and exponential 
infall rate. According to the assumption that there 
is lGyr delay for the halo formation, all of the gas 
infall throughout this paper begin at t = lGyr. 
It shows that, comparing with the star forma- 
tion rate predicted by Gaussian-infall model, the 
exponential-infall model predicts higher SFR at 
early stage and lower SFR later. This point di- 
rectly results in the differences of the mass-SFR re- 
lations at intermediate redshifts predicted by these 
two models (see following discussions). 

In order to compare with the observations, we 
plot the predictions of the exponential-infall model 
as dash lines in Fig. 4 and 5. It can be seen 
that, although both the solid and dash lines can 
be roughly consistent with the local observations, 
the exponential-infall model does predict lower 
SFR/M* than observations at intermediate red- 
shifts. Since the stellar population of local galaxy 
is resulted along its whole star formation history, 
local observations do not have the ability to distin- 
guish whether the model is viable at different time 
steps. Fortunately, the observed data at interme- 
diate redshifts are now available and can provide 
tighter constraints on the model of galactic evolu- 
tion. In other words, our results suggest that the 
Gaussian infall rate may be more reasonable to de- 
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scribe the gas cooling process than the exponential 
infall rate. 

In fact, by fitting synthetic energy distribution 
to the multi-band photometry of a sample of ac- 
tively star- forming galaxies at z ~ 2, Maraston et 
al. (2010) concluded that adopting an exponen- 
tially increasing SFRs and a high starting rcdshift 
appear to provide more reasonable results than 
that of an exponentially declining SFR in the case 
of star forming galaxies at redshift beyond ~ 1. 
Unfortunately, both their model and ours are phe- 
nomenological models. It is still beyond the scope 
of this paper to answer why the efficiency of gas 
cooling process or/and star formation process may 
be depressed at early stage and why the infall-peak 
time of low-mass galaxy is later than high-mass 
systems. 

Indeed, the semi-analytical models and stud- 
ies using hydrodynamical simulations may provide 
more physical and reasonable explanation on how 
gas accreted onto the galaxies than phenomeno- 
logical models. However, there is still a long way 
to go to fully understand the gas cooling processes. 
For example, the semi-analytical models assumes 
that gas within galaxy halos initially reaches the 
temperature of the virial halo, while some studies 
using high-resolution cosmological hydrodynami- 
cal simulations found that the majority of baryons 
may never reach the virial temperature as they ac- 
creted onto galaxies (Keres et al. 2005; Dekel & 
Birboim 2006). Brooks et al. (2009) found that 
low-mass galaxies are instead dominated by ac- 
cretion of gas that stays well below the viral tem- 
perature and even for galaxies at higher masses, 
cold flows dominate the growth of the disk at 
early times. All these findings called the general 
virial temperature assumption into question, but 
may give more information to improve the semi- 
analytical models and help us to understand more 
clearly on how the gas gets into the galaxies. 

Finally, comparing our model to that of BP00, 
the main difference is that BP00 adopted the expo- 
nential infall rate while we adopt the Gaussian in- 
fall rate. BP00 also notice the fact that, if the ex- 
ponential declining infall rate is assumed, the pre- 
dicted stellar population will be too old to repro- 
duce the observed blue colors of low-mass galaxies 
even if the infall time-scale tends to be as large as 
r —> oo. They assumed a negative infall time-scale 
for the most extended and less massive galaxies in 



their model, i.e. the infall rate increases with time. 
However, according to our discussions above, we 
guess that the model of BP00 may predict lower 
SFR/M* at intermediate redshift than the ob- 
served data although their model predictions can 
be in good agreement with the local observations. 

5. Summary 

In this paper, we present a phenomenological 
model for the formation and evolution of disk 
galaxies. It is assumed that the disk origins and 
grows by gas infall. We adopt the Gaussian form 
of gas infall rate with one free parameter, infall- 
peak time t p . We also include the contribution of 
gas outflow induced by supernova feedback. The 
Kennicutt star formation law is adopted to de- 
scribe how much cold gas turns into stellar mass. 
We calculate the evolution of disk galaxies and 
compare model predictions with observations. We 
find that: 

1. Model predictions are very sensitive to the 
infall-peak time t p in that the model adopt- 
ing later t p results in blue colors, low metal- 
licity, high specific star formation rate and 
high gas fraction, while gas outflow rate 
mainly influences the gas-phase metallicity 
and star formation efficiency mainly influ- 
ences the gas fraction. 

2. We construct a mass-dependent model by as- 
suming that low mass galaxy has later t p and 
larger b out than massive systems. It is shown 
that our mode predictions are in fairly agree 
with not only the local observations, but also 
the observed mass-SFR relations at interme- 
diate rcdshifts. 

3. It shows that, comparing with the star for- 
mation rate predicted by the Gaussian-infall 
model, the exponential-infall model predicts 
higher star formation rate at early stage and 
lower star formation rate later. Our results 
suggest that the Gaussian infall rate may be 
more reasonable to describe the gas cooling 
process than the exponential infall rate, es- 
pecially for low-mass systems. 

We emphasize that, here we present a phe- 
nomenological model. Through comparison be- 
tween model predictions and observations, our 
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main aim is to investigate main properties of the 
formation and evolution of disk galaxy, especially 
the trend as a function of galactic stellar mass. 
Unfortunately, there is still a long way to go to 
fully understand why low-mass galaxy may have 
later infall-peak time than high-mass systems and 
why the Gaussian infall rate may be more prefered 
than the exponential infall rate. 
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